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I.  INTRODUCTION  ] 

i 

i 

t 

The  assessment  of  delivery  accuracy  for  a weapon  is  usually 
initiated  through  test  firings.  The  results  of  these  test  firings 
are  miss  distances  of  the  rounds  from  the  target  center  in  the  j 
x (cross  range)  and  y (down  range)  directions,  usually  denoted  Ijy 

(x  , y^}  . These  n pairs  of  miss  distances  are  then  used  to  j 

i=l  ; 

compute  an  estimate  of  a parameter  or  parameters  which  described  the 
accuracy  of  the  weapon.  The  parameter  most  commonly  estimated  is 
the  CEP  (Circular  Error  Probable),  the  radius  of  a circle  withiti 
which  50%  of  the  future  rounds  will  fall.  This  estimate  provides 
a single  number  which  is  customarily  used  to  describe  the  accuracy 
of  the  weapon  under  study.  : 

The  procedure  described  above  can  be  improved  by  incorporating 
the  concept  of  statistical  tolerance  limits.  This  topic  is  addressed 
in  Reference  [8]  but  will  be  briefly  restated  here  for  completeness. 
Consider  the  case  where  the  miss  distances  from  the  target  centcir 
follow  a circular  normal  distribution.  This  distribution  is  giVen  by 


x -(xV)/(2o2) 
f(x,y)  a 2 e 

2nrT 


(1) 


where  n is  the  common  standard  deviation  of  the  miss  distances  ! 
(common  to  both  the  x and  y directions).  Under  this  distributional 
assumption,  the  relationship  between  CEP  and  ? is  well  known  to  be 

CEP  ■=  1.1774  o . ! (2) 


If  one  know  the  value  of  the  population  parameter  then  CEP  above, 
defined  in  terms  of  .* , provides  the  radius  of  a circle  which  contains 
uxactly  SQ%  of  the  bivariate  probability  and,  hence,  will  contain  50% 
of  the  future  rounds.  However,  ? is  usually  not  known  and  must  bo 
estimated  from  the  test  data.  One  such  estimate  based  on  the  theory 
of  maximum  likelihood  is  \ 


" » » . 5 ! 

‘ * { r <*:  + y:>/2ni  <?> 

i-i 


i 


where  the  cii*cumflex  (or  hat)  above  <j  is  used  to  signify  an  estimate 
of  the  true  population  parameter  <7.  Hence,  an  estimate  of  CEP  is 
obtained  by  replacing  a in  (2)  with  $ which  yields 


CEP  » 1.1774  S . 


(4) 


(This  estimate  is  slightly  biased  but  it  is  rarely  corrected  since 
the  bias  is  small  and  the  correction  factor  Ls  complicated.)  CEP 
above  is  only  a point  estimate  of  the  true  CEP,  and,  as  such,  it 
yields  no  information  regarding  the  probability  that  the  circle  it 
defines  contains  50%  or  more  of  the  population.  This  probability 
has  been  evaluated  in  Reference  [8]  for  the  circular  case  under 
discussion.  The  results  show  that  it  does  not  exceed  .50  for  any 
finite  sample  size.  This^means  that  one  has  at  most  50%  confidence 
that  a circle  of  radius  CEP  will  encompass  50%  or  more  of  the  future 
impact  points  from  the  weapon  under  test. 

To  increase  this  probability  or  confidence  to  a more  acceptable 
level,  it  is  clear  that  one  must  increase  the  multiplying  factor  in 
(4)  above  the  customary  1.1774.  Tables  of  these  multiplying  factors 
are  provided  in  Reference  [8]  for  the  case  of  circular  normally 
distributed  miss  distances.  These  factors  enable  the  analyst  to 
construct  circles  which  encompass  at  least  100P%  of  the  population 
(future  rounds)  with  100v%  confidence  for  a wide  variety  of  realistic 
values  of  n,  P,  and  v.  Circles  of  this  type  are  referred  to  in 
statistical  literature  as  tolerance  circles,  tolerance  regions  or 
tolerance  limits.  In  general,  (P,v)  tolerance  limits  are  functions 
of  Che  sample  data  and  define  a region  which  contains  at  least  10QP% 
of  the  population  with  lQ0v%  confidence. 

On  Che  basis  of  the  above,  one  can  utilise  developed  procedures 
So  construct  tolerance  circles  provided  Che  miss  distances  follow  • 
circular  normal  distribution,  lie  can  use  the  multiplying  factor 
which  will  provide  any  level  of  prescribed  confidence  not  only  for 
501  of  the  population  but  for  any  desired  percentage.  The  problem 
hero  is  that  miss  distance  distributions  are  not  always  circular 
normal.  In  fact,  mo  tv  often  titan  not,  they  tend  to  follow  an 
elliptical  distribution.  The  purpose  of  this  report  is  to  provide 
a procedure  for  constructive*  tolerance  cirelos  for  the  case  where 
the  miss  distance  distribution  follows  an  unear related  bivariate 
normal  distribution  with  unequal  variances  (hereinafter  referred 
to  as  elliptical  normal).  The  procedure  to  be  developed  in  the 
ensuing  sections  is  not  exact  in  the  sense  of  providing  circles 
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which  contain  at  least  100P7.  of  the  population  with  exactly  100y% 
confidence.  However,  the  confidence  afforded  by  these  approximate 
tolerance  circles  is  very  close  to  100y%.  The  closeness  has  been 
examined  through  a computer  simulation,  the  results  of  which  are 
provided  in  a later  section. 


II.  TOLERANCE  CIRCLES:  a and  0 KNOWN 

X Y 

It  would  be  instructive  at  this  point  to  examine  the  problem 

of  finding  the  radius  of  the  100P%  circle  when  the  miss  distances 

are  elliptical  normal  and  the  population  parameters  are  assumed 

known.  Such  circles  can  be  thought  of  as  tolerance  circles  when 

the  sample  size  n is  infinite  and  die  confidence  coefficient 

y - 1.00.  This  development  will  pave  the  way  and  ease  the  burden 

in  the  development  of  tolerance  circles  for  the  more  practical 

case  where  cv  and  n are  assumed  unknown. 

X Y 

We  shall  first  briefly  review  the  circular  case.  Let  the 
bivariate  random  variable  (X,Y)  represent  the  miss  distance  of  a 
round  in  the  x and  y directions,  respectively.  If  this  miss 
distance  follows  a circular  normal  distribution,  the  probability 
density  of  (X,Y)  is  given  by  equation  (1).  From  equation  (1),  it 

is  easy  to  show  that  the  probability  density  of  the  radial  miss 

x, 

2 2" 

distance  R » (X  + Y ) is 


%<t: 


*r2/2„2 


v > 0 


(5) 


where  the  subscript  C denotes  circular.  Furthermore . one  can 
obtain  the  cumulative  distribution  function  of  R in  closed  form, 
i.e., 


ec(r)  * i*rob(R  --  r)  » I ^<c)  dc  * i - a-*  ^ . <&) 

0 


■To  obtain  cho  well  known  rotation  between  -*  sm.*  tgp,  one  merely 
find*  the  value  of  t Ter  which  CL(r)  * .5.  Solving 


one  finds  r = CEP  = 1.1774?  which  was  given  in  equation  (2).  To 
find  t*.is  radius  of  til c general  1001*%  circle,  C^,  one  solves 


P = 1 


for  r wlsich  yield* 


v 


= a f-  2 In  (l-£>)}‘. 


(7) 


The  development  is  quite  easy  for  the  circular  case  since  only  one 
population  parameter,  t,  is  involved  and  the  cumulative  distribution 
function  of  R can  be  expressed  in  closed  form.  Let  us  now  turn  to 
^he  elliptical  case. 

If  the  miss  distance  follows  an  elliptical  normal  distribution, 
the  density  of  (X,Y)  is  given  by 


-f  (*/<?,.)“  + Cy/s?v)2}/2 

t'(x.y)  * (1/2 * -«s  y)  e * ' (S) 


where  «„  and  are  the  raisa  distance  standard  deviations  in  the 
s and  y directions , respectively,  The  probability  density  &i  the 
radial  e»i*s  distance  R is  wm  complicated  in  this  case  and  is  given 

fey 


ggCf)  “•< y)  e ly(brta),  r ? 0 (9) 


where  the  subscript  £ denotes  elliptical  anJ  where 


* 2 
V -V 


b 


■< 


S 


ami  I lx)  i9  the  Modified  Bessel  Function  of  the  First  Kind  end  Zero 
Order,  i.e., 

>0w%Kx“‘8«- 

0 

The  derivation  of  the  density  of  R in  <9)  is  obtained  by  applying 
the  usual  transformation  from  Cartesian  to  polar  coordinates  to 
(8).  The  complete  derivation  is  provided  by  Chew  and  Boyce  in 
Reference  [l]. 

The  cumulative  distribution  function  of  R is  defined  by 

r 

GE(r)  = Prob(R  <;  r)  = J S£(t)  dt.  (10) 

0 

However,  it  is  clear  that  Gt>(r)  cannot  be  expressed  in  closed  form 
since  the  integrand  contains  an  integral  which  cannot  be  expressed  in 
closed  form.  Hence,  the  radius  of  the  100P7,  circle  for  the  elliptical 
case  cannot  Ik?  expressed  by  a simple  formula  as  it  was  for  the  circular 
case  (equation  (?)).  thw  muse  solve 


l*  « a.(r> 


HD 


for  r - t'  which  can  only  be  aecowplUhed  by  nvaawrical  methods. 

The  direct  solution  of  equation  (U)  (solution  of  P for  fixed 
values  of  r,  ami  *y)  can  be  obtained  hy  using  ao  efficient 
algorithm  developed  by  hlheoat©  and  iarnaoin  clutch  is  dotusented 
Reference  i ■ . The  r.o Kjt  {*•'.»  (solution  of  r for  fix«d 

vy|«e*»  of  p,  .jj,  -y>  can  bo  obtained  item  table*  prepared  by 
Pi8©«««*  and  dWnasltt  and  set  forth  in  the  sasw  document.  The 
(mint  to  be  -iado  f*  chat  to  find  ih  - radius,  of  the  tWW  circle  for 
the  elliptical  c4*o,  one  rs*«*ft  to  table*  ©r  engage  in  ©xtetutive 
pueterltai  effort*. 

thy*  far,  no.  *e»tl©n  ha*  ho«n  about  the  (j|p  tor  the 
elliptical  ■>' .4*0.  The  tva>en  ss  that  the  c£r  per  1*  associated 
rn.lv  «UH  the  circular  esse,  however.  is  certainly  a ©title 

centered  sf  the  target  center  ©tilth  contain*  %t?Z  of  the  bivariate 
probability.  The  radius  of  the  circle  is  jw*t  a Special  case  of  C& 


above  with  P -•  .50.  To  distinguish  it  from  the  CEP  for  the  circular 
case,  the  radius  of  the  50%  circle  for  the  elliptical  case  is  some- 
times referred  to  as  an  "equivalent  CEP"  (ECEP) , The  reason  for 
distinguishing  between  the  two  is  to  avoid  possible  misuse  of  the  CEP 
for  the  elliptical  case.  For  example,  if  one  were  given  a CEP  with 
no  mention  that  it  pertained  to  the  elliptical  case,  he  may  be  tempted 
to  divide  it  by  1.1774  and  use  the  resulting  a in  equation  (6)  to 
obtain  "hit  probabilities."  Such  probabilities  are  erroneous  for 
values  of  r # ECEP  unless  the  population  distribution  is  circular 
normal. 

To  avoid  a table  "look  up"  for  the  ECEP.  several  approximations 
have  become  fashionable  over  the  years.  The  three  most  common  are 
given  below: 


ECEP  « 1.1774  (ax  ay)^ 

(12) 

ECEP  ~ 1.1774  [ (ux  + ay)/2] 

(13) 

ECEP  ^ 1.1774  [(ax  + aY)/2]  . 

(14) 

These  approximations  are  discussed  by  Groves  in  Reference  [4]  where 
it  is  shown  that  unless  ^x  is  very  close  to  o on^y  (13)  has  merit 
and  then  only  if  the  ratio  of  the  larger  to  tne  smaller  of  these  two 
standard  deviations  is  less  than  about  five.  A fourth  approximation 
will  be  given  later  in  this  report  as  a limiting  form  of  the  expression 
to  be  derived  for  tolerance  circles.  This  fourth  approximation  is  not 
restricted  to  the  radius  of  the  50%  circle  but  has  general  application 
to  the  10QP%  circle.  It  is  found  to  be  more  accurate  than  any  of  the 
above  three  for  approximating  the  radius  of  the  50%  circle. 

Throughout  this  section,  it  has  been  assumed  that  a x and  are 
known  a priori.  This  is  rarely  the  case  in  that  and  oy  are  usually 
estimated  from  test  data.  One  could  use  these  estimates  in  place  of 
ojj  end  oy  in  equation  (11)  to  obtain  point  estimate  of  the  ECEP  or  tho 
more  general  Cp.  However,  point  estimation  for  the  elliptical  case 
yields  no  more  information  than  for  the  circular  case.  Therefore,  lot 
us  turn  our  attention  to  the  problem  of  deriving  tolerance  limits  for 
tho  elliptical  case  when  the  parameters  are  estimated  from  tost  data. 
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III.  APPROXIMATE  TOLERANCE  Cl ROLLS : ax  AND  ay  UNKNOWN 

Suppose  n rounds  are  fired  at  a target  and  yield  miss  distances 

fx  , y. }n  . Under  the  assumption  that  the  miss  distances  follow  an 
1 1 i=l 

elliptical  normal  distribution,  the  probability  density  for  the 
population  of  miss  distances  is  given  by  equation  (8).  The  standard 
deviations  and  are  unknown,  but  their  maximum  likelihood 
estimates  are 


= f E xf/n] 
'i-1  1 


<?Y  " { £ y-/“} 

* i=i  1 


To  construct  a tolerance  circle  in  accordance  with  the  definition  on 

page  3,  we  need  to  find  a function  of  thr  sample  data,  say  L(x^, 

x > y, , . ..,  y ) which  satisfies 
n 1’  ■'n 


J I zri 


-[ (x/ox)2  + (y/aY)2]/2 


2^  2^t  2 

X 4y  sL 


°X  aY 


dxdy  2 P>=  y (15) 


where  the  arguments  of  L have  been  deleted.  In  the  above,  L is  the 
radius  of  a circle  which  contains  at  least  1Q0P%  of  the  population 
with  100y%  confidence.  It  is  easy  to  show  that  ax  and  <jy  above  are 
sufficient  for  their  respective  parameters.  Hence,  the  radius  L of 
the  tolerance  circle  can  be  expressed  ns  a function  of  these  two 
estimates  vice  all  2n  observations.  That  is,  we  can  express  L in 
functional  form  as  L(,Tj.,  oy)  vicu  L(Xj,  x^,  y^,  . ..,  y^). 

Equation  (15)  is  a general  expression  of  the  problem.  However, 
the  problem  can  be  expressed  more  simply  since  a tolerance  circle  for 
this  case  is  equivalent  to  an  upper  tolerance  bound  for  the  distribution 
of  the  radial  error  K as  given  by  equation  (9).  Hence,  tho  problem  can 
bo  restated  as  one  oi  finding  L(<\. , &.,)  which  satisfies 

« I 


8 


Prob  J (r/axaY)  e"  ar  IQ(br2)  dr  ;>  P J = Y 
0 


(16) 


where  a and  b are  the  functions  of  rrx  and  a-v  defined  earlier.  The 
function  L(a^,  <Ty)  above  could  be  found  if  (a)  the  integral  in  (16) 
could  be  expressed  in  closed  form  (as  in  the  circular  case)  or  (b) 
a suitable  transformation  could  be  found  to  express  the  integrand 
free  of  the  unknown  population  parameters  Oy  and  ay.  Condition  (a) 
clearly  does  not  hold  in  this  case,  and  with  respect  to  condition 
(b),  a suitable  transformation  has  not  been  found.  Hence,  an 
approximate  solution  will  be  sought. 

The  integrand  in  £16)  is  the  distribution  of  the  radial  miss 
distance  R = (X^  + Y^)'2  under  the  assumption  of  elliptical  normality 
on  the  bivariate  random  variable  (X,Y).  Consider  now  an  approximation 
to  the  distribution  of  a function  of  R.  Under  our  distributional 

•y  •y 

assumption,  X ~ N(0,  ox),  Y ~ N(0,  o“),  and  X and  Y are  independent 
where  signifies  "is  distributed  according  to"  and  N(p,,a^)  signifies 
a normal  distribution  with  mean  p and  variance  a^.  It  is  well  known 
that  the  square  of  a standard  normal  variable  (N(0,1))  is  a chi-square 
variable  and  that  independent  chi-square  variables  are  additive  where  the 
density  of  a chi-square  variable  with  v degrees  of  freedom  (y2)  is  given 
by 

h(w)  » t x w(v/2>  1 e W^2,  w 0 . (17) 

r(v/2)  2vU 


2 2 2 

Since  X and  Y are  not  standard  normal  variables,  R = X +Y  is  uot 
a chi-square  variable.  However,  the  distribution  of 


U 


o 


(18) 


can  be  approximated  with  a chi-square  density  with  v degrees  of 
freedom  where 


2 2 ' 
(r>x  + Oy) 

“ . 4 

’X  + 


(19) 


This  approximation  is  a special  case  of  a more  general  approximation 
due  to  Satterthwaite  (see  Reference  [6]).  Tha  general  form  is  as 
follows: 


Let  Vi  X^/a?  ~x^  > 1 = 1 k 


and  let  X^,  X^  be  independent. 


Then  U = 


k 

v E 
i=l 
k 

E a 
i=l 


a. 

1 


i 


2 


(20) 


approximately  where  the  a^  are  constant  and 


I k 2 

E ai  ci 
i=l 


k 

v 

i=l 


(ai  °Vv±) 


For  our  application  above,  k 2,  a^  = a..^  * 1,  and  55  ^2  a 1« 

The  approximate  distribution  of  U in  (18)  could  be  used  to  find 
the  approximate  distribution  of  R,  However,  one  can  restate  the 
problem  in  terras  of  H via 


A 

Probfj*  h(u)du  ;»  P]  « v (21) 

0 
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2 2 2 

where  A * \>L  / (a^  + cry)  and  h(u)  is  the  chi-square  density  given  in 
(17)  with  degrees  of  freedom  v given  by  (19).  The  inequality  in  (21) 
holds  if  and  only  if  A > p where  y2  p is  the  100P%  poin*-  for  a 
chi-square  variable  with  v degrees  of  rreedora.  That  is,  if  U ~ ^2, 
then  y2  p j[s  such  that  v 


Prob  (U  < ^p)  = P . 


Hence,  (21)  can  be  restated  as 


(22) 


or  equivalently, 


(23) 


To  complete  the  solution,  we  need  to  obtain  the  upper  100y%  confidence 
bound  for  ay  + cry  • We  then  equate  this  bound  to  VL 2/v2  D and  solve 
for  L.  v,r 

2 2 

To  obtain  a confidence  bound  for  ov  + oy,  we  again  resort  to 
Satterthwaite's  approximation  (20).  Under  the  initial  normality 
assumption, 


„*2/  2 2 . .2,  2 2 

aX/aX  ~ \\  and  \ * 


11 


Using  (20)  with  k - 2.  ai  . s,  ■ 1,  and  Vl  . di!trlbutl(m 


v '(<?x  + <7y ) 

rnr- 

°X  + °Y 


offreedorwhiS"^  * Chi'SCIUare  distribution  with  v'  degrees 


/ 2 , 2 * 

,/  .i2L LfxL 

4.  4 “ nv 

ax/n  + CTy/n 


Using  this  approximation,  the  upper  100Y%  confidence  bound  for 
ax  + cjy  is 


We  now  set 


nv(nx  + ay) 
2s— 

*riv,l-y 


vli  _ nv($x + 
2 2 


and  solve  for  L.  We  find 
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/ 


(24) 


L (ax,  $Y) 


2 2 ^ 
<ax  + aY) 


2 2 

The  parameter  v in  (24)  is  a function  of  the  parameters  and  cry 
which  are  unknown.  Hence,  v must  be  estimated  from  the  data  by 
replacing  with  and  with  aY  in  formula  (19)  for  v.  One 
notes  that  v is  not  restricted  to  the  integers  in  (24)  so  that  to 
apply  this  formula  one  needs  a table  of  chi-square  percentage  points 
for  fractional  degrees  of  freedom.  A table  of  this  type  has  recently 
been  prepared  by  DiDonato  and  Hageman  (Reference  [2])  and  will  be 
utilized  in  the  example  which  follows. 

Suppose  15  rounds  are  fired  at  a target  and  provide  the  miss 
distances  shown  in  Table  1.  Estimates  of  oj-  and  oY  obtained  from 
the  data  are  <?x  = 85.11  and  <$Y  = 20.55.  Let  us  use  these  values 
to  construct  (P,v)  tolerance  circles  for  P = .50,  y = .90  (the 
circle  which  contains  50%  of  the  population  of  future  rounds  with 
90%  confidence)  and  for  P = .50,  v - .95  (the  circle  which  contains 
50%  of  the  population  of  future  rounds  with  95%  confidence).  For 
both  circles,  the  value  of 

($1  +nY)‘  = [(85.il)2  + (20.55) V = 87.56. 


we  now  need  to  calculate  v which,  according  to  formula  (19)  is 


a ! (85. nr  + ao.ssrV 

(85. 11)*  + (2Q.55)4 


1.116  . 


and  itv  “ 15v  a 16.74.  The  value  of  the  507'  (value  of  P)  point  is 
needed  for  a chi-square  with  1.116  degrees  of  freedom  as  well  as 
values  of  the  5t  and  10%  (values  of  l*v)  points  for  a chi-square 
with  16.74  degrees  of  freedom,  Reference  [V  cited  above  provides 
tabular  values  of  chi-square  percentage  points  for  *•  * .1C.1) 
13(.12)2&(*25)42.5(. 5)512.  The  points  needed  for  our  degrees  of 
freedom  can  be  obtained  using  those  table*  with  linear  Interpolation, 
or  they  can  be  obtained  exactly  by  using  the  program  documented  in 
the  report.  Linear  interpolation  of  tabular  value*  provide* 
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TABLE  1 


Hypothetical  Miss  Distances  For  Sample  Problem 


x 

11.71 

51.94 

- 87.60 
31.11 

-104.12 

23.95 

- 15.10 
96,24 

157.33 

- 74.22 
68.04 

- 90.03 
-110.75 

95.01 

111.74 


« 85.lt 


—1 

- 3.40 
-19.68 
-20.80 
12.35 
15.68 
1.53 
16.60 
-12.55 
-47.80 
3.43 
-14.78 
-24.70 
8.05 
-21.05 
32.15 


- 20.55 


.5564 


X1 .116,  .50 


x16.74,  .10 


9.8836 


2 

x16.74,  .05 


8.4863  . 


The  exact  values,  obtained  by  running  the  program,  are  .5564,  9,8822, 
and  8.4844,  respectively.  The  values  obtained  by  linear  interpolation 
are  the  values  which  would  usually  be  available  to  the  analyst.  They 
arc  sufficiently  accurate  and  will  be  used  to  complete  our  sample 
problem. 

For  (P,v)  **  (.50,  .90),  we  find 


- x -CL5IG55k-aV" 

Y;  { 9.8836  j 


(87.56)  = 80.46, 


and  for  (l»»v)  * 1.50,  .95),  we  lino 


3C 


{ 


(ISM. 5564)  \ * 
8.4863"  j 


(87.56) 


86.83  . 


The  point  estimate  m the  FtdvP  for  this  case  is  61.6?,  It  is 
obtained  by  outre tny  the  tables  In  Reference  f ')!  with  e » 20.55*85.11 
* .2415  and  l'  « .51)  to  obtain  .7217  which  when  multiplied  by  ^ 
yields  tlu-  above  result,  it  has  an  associated  confidence  of  teas 
than  50&,  i.e.,  *<r  are  less  than  S0£  ewfiJent  that  a circle  of 
radius  61.42  contains,  at  lease  SOI  of  the  population  of  impact*. 

The  data  points  for  this  oxa-uple  and  all  three  circles  are 
shown  in  Figure  l.  All  three  circles  can  be  thought  of  as  estimate* 
of  the  "50t  circle**  hut  with  different  levels  of  confidence,  One 
notes  in  this  example  and  In  general  that  the  cenfidenee  increase* 
a*  the  ci.cle  radius  increase*. 


IS 


f 


* 


FIGURE 


IV.  ACCURACY  EVALUATION 


Satterthwaito's  approximation  (20)  was  used  twicf'  during  the 
formulation  of  the  (!’,v)  tolerance  circle  given  by  equation  (24). 
Therefore,  it  is  an  approximate  tolerance  circle,  and  one  should 
know  the  worth  of  this  approximation  in  order  to  use  it  with 
"confidence."  The  purpose  of  this  section  is  to  evaluate  the  accuracy 
of  the  approximation;  several  approaches  will  be  used  in  this  endeavor. 
Let  us  first  define  0 •*  c = ry/ 'x  ' Then  the  first  approach  will  be 
to  examine  L(^v»  in  (?A)  for  the  cases  where  c = 1 and  c « 0.  This 
will  correspond,  respectively,  to  the  circular  case  (case  where 
crx  ==  rfy)  and  the  univariate  case  (case  where  rr^  differs  so  much  from 
("Y  that  the  entire  population  is  concentrated  on  a line).  The  second 
approach  will  he  to  examine  the  confidence  afforded  by  the  approximation 
for  small  sample  sizes  n using  a Monte  Carlo  simulation.  The  last 
approach  will  be  to  examine  the  limit  of  L(ffx,  try)  as  n increases 
without  bound.  The  results  of  these  examinations  should  be  sufficient 
to  pass  judgment  on  the  worth  of  equation  (24). 

If  c = l,  then  *x  - -y  a t,  say,  and  rx  and  in  (24)  are  both 
estimating  the  common* variance  so  that  oZ  + becomes  2S^,  Also, 
v,  in  equation  (IV)  becomes  * ' 


•> 


Hence,  If  e * 1.  Ureses 


it 


Equation  (25)  coincides  uxactly  with  the  exact  solution  for  tolerance 
circlos  for  the  circular  casu  provided  by  Thomas  and  Crlgler  in 
Reference  [7].  (Reference  [7]  is  a refinement  of  the  work  in 
Reference  [8]  for  publication  in  the  open  literature.)  Therefore, 
if  ■ fTy,  the  approximate  solution  in  (24)  corresponds  to  the 
exact  solution  for  the  circular  case. 


-2 


To  examine  c ■ 0,  we  shall  simply  set  <jy  " 0.  Here  we  find  + 
and  v in  equation  (19)  equals  one  soYthat  L(a^*  <5y)  becomes 


L«x> 


(26) 


Equation  (26)  is  the  exact  expression  for  the  upper  tolerance  bound 
for  a univariate  normal  distribution  when  the  mean  is  known  a priori 
to  be  zero.  (See  the  development  of  the  univariate  case  by  Hald,  for 
example,  in  Reference  [5].)  Therefore,  the  approximate  solution  in 
(24)  becomes  the  exact  solution  in  the  degenerate  case  where  one  of 
the  standard  deviations  is  zero. 

The  results  of  the  last  two  paragraphs  indicate  that  the 
approximation  provides  what  is  needed  for  values  of  c close  to 
zero  or  one.  For  mid-range  values  of  c,  we  shall  have  to  resort 
to  Monte  Carlo  simulation.  A single  replicate  of  the  simulation 
process  used  in  this  study  is  described  below: 

(i)  For  given  values  of  P,  y,  ct^»  and  ay, 
a sample  of  size  n was  generated  from 
an  elliptical  normal  distribution  using 
Monte  Carlo  sampling  techniques. 

(ii)  Using  the  data  points  generated  in  (i), 

6j(,  cTy>  and  v were  computed. 

(iii)  Using  the  value  of  v computed  in  (ii) 
and  the  specified  values  of  P and  y, 
the  required  fractional  chi-square 
percentage  points  were  computed  and 
L(a^,  6y)  formed. 
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(iv)  The  Integral  I of  the  elliptical  normal 
distribution  with  parameters  ax  and  oy 
was  computed  over  a circle  of  radius 
b(nx,  fly). 


The  above  process  was  repeated  N times,  and  the  confidence  coefficient 
was  estimated  by  the  ratio  of  the  number  of  replicates  In  which  1 > P 
to  the  total  number  of  replicates  N.  If  the  equation  for  L were  exact, 
then  y would  lie  within  sampling  variation  of  y (the  estimate  of  the 
true  confidence)  for  all  values  of  P,  y,  c,  and  n.  Since  L in  (24) 
is  an  approximation,  the  departure  of  y from  y (outside  sampling 
variation)  is  a measure  of  the  worth  of  the  approximation. 

The  selected  values  of  the  input ' parameters  P,  y,  n,  and  c for 
the  simulation  wore  as  follows: 


P:  .50,  .90 
y:  .90,  .95 
n:  5,  10,  20 

c:  0,  .05,  .10,  .20,  .25,  .33,  .50,  .57,  .67,  .80,  1.00. 

Values  of  n were  restricted  to  small  samples  in  order  to  avoid 
excessive  computer  running  time.  The  number  of  replicates  N for 
each  combination  of  parametric  values  was  set  at  10,000.  This  value 
insures  (with  957.  confidence)  that  the  true  confidence  (not  necessarily 
y)  lies  within  .01  of  v*  The  results  are  shown  in  Figures  2-5  vhere 
y is  plotted  as  a function  of  c in  multiple  curves  representing  the 
three  different  values  of  n.  Each  figure  represents  a different  (P,y) 
combination.  For  each  combination,  horizontal  lines  have  been  drawn 
at  v ± .01.  Values  of  y within  the  region  enclosed  by  these  lines 
can  be  considered  within  sampling  variation  of  y.  Values  outside  the 
region  represent  values  of  y outside  sampling  variation  and  indicate 
departures  of  the  approximation  from  y.  One  notes  that  in  all  four 
cases  y lies  in  the  region  of  sampling  variation  for  c close  to  zero. 
This  merely  confirms  the  results  of  the  earlier  discussions  on  the 
behavior  of  the  approximation  when  ov  and  oy  are  substantially 
different.  One  does  not  observe  a similar  behavior  for  all  values  of 
P,  y,  and  n for  the  case  where  c * 1.  He  does  observe,  however,  that 
for  the  larger  value  of  n(n  ■ 20),  y lies  either  within  or  very  close 
to  the  region  of  sampling  variation.  Departures  also  exist  for  the 
mid-range  values  of  c (most  notably  for  values  of  c close  to  .25)  with 
the  departure  being  less  sovere  for  the  larger  n than  for  the  smaller 
one. 
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06' 


20 


0 6 


/<r 


On  the  basis  of  Figures  2-5  and  the  above  comments,  one  can 
conclude  that  the  approximation  is  very  good  for  values  of  c close 
to  zero  regardless  of  the  values  of  1’,  y,  and  n.  For  values  of  c 
close  to  one,  the  approximation  is  also  very  good  if  n is  not  too 
small.  For  other  values  of  c,  there  are  departures  in  the  confidence 
afforded  by  the  approximation  which  appear  to  increase  as  P increases 
and  diminish  as  n increases.  We  now  need  to  insure  that  the 
approximation  continues  to  improve  as  n increases  beyond  the  maximum 
value  of  20  used  in  the  simulation.  Hence,  we  shall  consider  the  limit 
of  L(<Tj.,  ay)  as  n -.  jo  . it  is  easy  to  show  that 


1 


The  right  hand  portion  of  the  limiting  expression  is  based  on  the 
fact  that  ^2  and  ~~  are  consistent  estimators  and,  hence,  approach 
and  respectively,  as  n - -r,  The  left  hand  portion  is  based 
on  the  fae*  that  as  the  degrees  of  freedom  f — the  distribution 
of  xjr/f  is  concentrated  at  unity.  Hence,  in  the  limit,  all  percentage 
points  of  a chi-square  random  variable  divided  by  its  degrees  of  free- 
dom are  unity.  If  approximation  (24)  continues  to  improve  as  n - 
chon  the  percent  of  the  population  encompassed  by  a circle  with  radius 
given  by  (27)  should  approach  P.  The  percent  of  the  population 
encompassed  by  circles  of  radius  (27)  was.  evaluated  for  e = Q.0(, 1)1.0 
and  p a .50,  ,40,  .95,  and  .09,  The  results  are  set  out  in  Table  2 
and  show  that  while  those-  percentages,  designated  l* *,  are  pot  exactly 
equal  Co  t\  they  ate  very  nearly  equal  to  P,  differing  only  in  the 
3rd  decimal  place,  the  meaning  here  is  that  Che  percent  of  the 
population  encompassed  by  a circle  oi  radius  b(o^,  *y)  -»  P # vice  P 
as  n ^ *.  This  pesos  no  real  problem  since  P*  is  close  to  i\ 

Hence,  one  can  rest  assured  that  for  all  practical  purposes, 

!.(?£,  o«)  in  (22*)  improves  as  n increases  without  bound. 

We  can  sumartro;  the  worth  of  approximation  (24)  with  the 
following  remarks; 

(a)  It  -v, -preaches  the  exact  value  as  e « 0 
and  appears  to  produce  excel loot  results 
"regardless  of  a,  P,  and  v for  values  of 

e cUy»?  to  s-arv. 


limit.  L(A  ,rt  ).'»( 

II  -1  .TO 


(a.,  + 


». 

2* 

ctv) 


(27) 


TABLE  2 


Evaluation  of 


■>2  2 
x“4y  <-(1  imit  L) 


-[  <x/ax)2+(y/ov)2]/2 
o A * dxdy 


c 

JL 

.0 

.50 

.5000 

.1 

.50 

.5011 

.2 

.50 

,5032 

.3 

.50 

, 5018 

.4 

.50 

,4981 

.5 

.50 

,496? 

.6 

.50 

. '+964 

.7 

.50 

.<+976 

.8 

.50 

. ! 9K9 

.9 

.50 

.<‘<997 

1.0 

.50 

. -<000 

.0  - 

.95 

. 95  00 

,1 

.95 

.9500 

.2 

.95 

.9501 

.3 

.95 

,9503 

.4 

.9506 

.5 

.95 

.9510 

.6 

.95 

.9513 

.7 

.35 

,**5U 

.8 

.95 

.9504 

.9 

.93 

.9502 

1.0 

.95 

.9500 

c 

JL 

P' 

.0 

.90 

.9000 

.1 

.90 

.9004 

.2 

.90 

.9014 

.3 

.90 

.9028 

.4 

.90 

.9042 

.5 

.90 

.9052 

.6 

.90 

.9049 

.7 

.90 

.9034 

.8 

.90 

.9017 

.9 

.90 

.9004 

1.0 

.90 

.9000 

.0 

.99 

.9900 

.1 

,99 

.9899 

.2 

. .99 

.9897 

.3 

.99 

.9893 

.4 

.99 

.9890 

.5 

.99 

.9889 

.6 

.99 

,9890 

.7 

.99 

.9893 

.8 

.9897 

.9 

.99 

.9899 

1.0 

.99 

.9900 

(b)  It  approaches  the  exact  value  as  c -*  1, 
but  the  results  do  depend  on  n. 

(c)  For  intermediate  values  of  c,  there  is 
some  departure  from  the  desired  confidence} 
but  even  for  a sample  size  as  small  as 

n = 5,  it  does  not  exceed  3%  and  is  much 
less  for  most  values  of  c. 

(d)  It  improves  for  all  values  of  the 
parameters  as  n increases. 

Hence,  the  approximation  appears  to  be  suitable  for  general  applications 
with  little  loss  in  accuracy  even  for  sample  sizes  as  small  as  n - 5. 


V.  MULTIDIMENSIONAL  EXTENSION 


Thus  fiir  our  attention  has  been  confined  to  the  two  dimensional 
case  where  miss  distances  arc  recorded  in  the  (x,y)  plane.  This 
could  be  the  ground  plane  or  perhaps  the  plane  normal  to  the  trajectory 
at  impact.  In  the  case  of  an  air  burst  weapon,  however,  the  miss 
distances  are  throe  dimensional,  and  for  other  applications,  the  number 
of  dimensions  could  exceed  three.  In  these  cases,  we  would  be 
concerned  with  a k dimensional  tolerance  sphere  vice  a tolerance 
circle.  Extending  the  work  of  the  previous  sections  to  tolerance 
k-spheres  for  the  multidimensional  case  is  straightforward  and  will 
be  outlined  below. 

LcL  the  multivariate  random  variable  (Xj,  ...,  X^)  represent  the 
miss  distance  of  a round  in  k directions.  If  this  miss  distance 
follows  an  uncorrelated  multivariate  normal  distribution,  its 
probability  density  is  given  by 


1 (*i| i •••,  ) 


(2n> 


k/2 


-%r  (x^/rf >“  + ...  + (Xj./(7k)23 
e (28) 


l ^2  *'  * 


where  is  the  miss  distance  standard  deviation  in  the  ich  direction. 
If  a sample  of  n rounds  is  fired,  the  n miss  distances 

...»  xkj’f*l  are  used  to  estimate  the  <s^.  The  method  of 
maximum  likelihood  yields 


i ” * i 

a i £ sT  Jn)  . i * l,  2, 

l*-‘  J J 


• * * * 


To  construct  a tolerance  k- sphere  we  need  to  find  the  radius  L as 
a function  of  the  sueh  that 
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where  the  region  of  integration  S is  the  interior  region  of  the 
k-sphere  xj  + ...  + xj*  = 1,2.  If  wo  let  r2  « X2  + ...  + X2  tlie 
application  of  (20)  provides  that  1 k 


U = 


CT1  + ...  + ak 


2 ~ \ approximately 


where 


2 

(cij  + 


4 + 


“ +^k> 


► H 

>•  + at< 


(30) 


One  can  now  restate  (29)  in  terms  of  the  distribute -a  of  U via 


A 


l* fob  1 J h(u)du  P 


(31) 


0 


where  h(u)  is  the  chi-square  density  with  degrees  of  freedom  v given 
by  (30)  and  A » vh“/{^  + ...  + ^),  Equation  (31)  holds  if  ^niy 

if  A > p*  Neneo,  (31)  can  be  rewritten  in  fora 


Prob 


vl‘ 


2 2 

^ ^ k 


er  equivalently 


Applying  (20)  again,  we  find  the  upper  100y%  confidence  bound  for 
Oj  + ...  + crj-  to  be  approximately 


.2  2 

nv  (oj  + 5k) 

! 


inequality  in  (32),  one 


(c*  + ...  + <?£)'  . (34) 


One  notes  that  this  general  expression  is  very  similar  to  the  two 
dimensional  expression  given  by  (24).  The  difference  lies  in  the 
formula  for  v and  the  additional  terms  in  the  right  hand  portion 
to  account  for  the  increased  dimensionality.  This  expression  for 
L is  the  radius  of  a k dimensional  sphere  which  contains  at 
least  1Q0P%  of  the  population  with  confidence  approximately  equal 
to  lOOvE.  The  worth  of  this  approximation  has  not  been  studied  for 
dimensionality  above  two. 


*nv,l-> 


Equating  (33)  to  die  right  hand  side  of  the 
obtains 


, r ~ v _ / %.P 

Hcr^,  ....  cTj. > - \ y* — 

\iv,l-y 


(33) 


VI.  ANOTHER  APPROXIMATION:  PARAMETERS  KNOWN 


j 

| 


It  was  shown  in  Nation  II  that  when  the  parameters  are  known 
a priori,  the  ECEP  or  the  radius  of  the  general  100P7.  circle  for  the 
elliptical  case  can  be  obtained  exactly  but  only  by  table  look-up  or 
extensive  calculations.  Several  approximations  were  shown  ((12),  (13), 
and  (14))  which  are  oftentimes  used  for  KCEP  calculations  to  avoid  the 
exact  procedure.  A fourth  approximation  which  can  be  applied  to  the 
general  caSo  of  the  lOOPft  circle  is  simply  the  limiting  fora  of 
L(r?x,  jy)  In  (24)  as  n — -o.  This  form  was  given  in  equation  (27)  in 
Section  IV  and  upon  minor  rearrangement  becomes 


1 J-,  t. 

2\,,P/vf  [(<TX  + aY)/2}2  (35) 


where  v is  given  by  (19). 

thie  notes  from  Table  2 that  this  approximation  is  exact  for  the 
degenerate  cases  where  c “ l*  Also  from  Table  2,  wc  see 

that  for  midr.utge  values  oi  c,  a circle  of  radius  (35)  contains  very 
nearly  1O0P?.  c£  the  population,  differing  from  P wily  in  the  3rd 
decimal  place.  However,  to  employ  approximation  (35),  we  need  values 
of  \~  for  fractional  v.  To  obtain  such  values,  we  have  to  resort 
to  tables  or  engage  *n  heavy  calculations.  These  are  the  same 
alternatives  we  faced  for  the  exact  solution.  Therefore,  one  could 
argue  that  it  the  approximation  involves  the  use  of  tables  of 
percentage  points  for  the  chi-square  with  fractional  degrees  of 
freedom  (Re£erom;e.  * ;!),  why  not  use  inverse  tables  of  the  circular 
coverage  function  (Reference  f 3})  to  obtain  the  exact  solution?  Tins 
Ik  valid  argument,  jnd  no  attempt  will  be  made  to  refute  it  except 
to  sav  that  many  Analysts  feel  more  comfortable  using  a multiplying 
l actor  for  so-v  function  of  and  xO  than  they  do  using  the  table 
look-up  for  the  esWt  solution,  furthermore , if  one  confines 
his  attention  tn'rlie  EgEP,  lie  can  readily  construct  a short  table  of 
such  factors  for  oayv  access . Table  3 below  is  a short  table  of  this 
type  where  the  percentage  points  have  been  taken  from  Reference  f 2] . 

To  Illustrate  the  use  of  Table  3 and  to  compare  approximation 
(IV)  with  the  exact  xelut ion  and  the  other  approximations,  consider 
the  lot  lowing  example;  If  »*y  **  15  ft*  and  **  30  ft. , 
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TABLE  3 


2 2 

Multiplying  Factors  of  [ + Cy)/2] 
Lor  ECEP  Approximation 


iL&joM 


1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 
1.9 
2.0 


.9538 

.9928 

1.0258 

1.0542 

1.0789 

1.1005 

1.1195 

1.1365 

1.1516 

1,1652 

1.1774 

2 , 2 
(BX  * 

B it-  ■■■■ 


find  the  BCKP  by  the  exact  method  and  compare  with  the  four 
approximate  solutions.  We  obtain  the  exact  solution  frees  tables 
in  Reference  13}.  The  first  three  approximations  are  straight  forward 
and  require  no  elaboration,  the  last  approximation,  (35),  involves 
the  use  of  Table  3.  For  this  ease 


u B K30)2  * qs?5i 

(30)4  + (IS)4 


1.4704. 
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Interpolating  between  v = 1.4  and  1.5  in  Table  3 yields  a multiplying 
factor  equal  to  1.0941.  The  results  are  set  out  below  for  comparison. 


f1y  = 15  ft.. 


.8704 


frx  = 30  ft. 


Exact  ECEP  « .8704  ? = 26.1120  ft..  Reference  f3] 

Approx.  ECEP  = 1.1774(21.2132)  = 24.9764  ft..  Equation  (12) 

Approx.  ECEP  = 1.1774(22.5000)  = 26.4915  ft..  Equation  (13) 

Approx.  ECEP  = 1.1774(23.7171)  = 27.9245  ft..  Equation  (14) 

Approx.  ECEP  = 1.0941(23.7171)  ^ 25.9489  ft..  Equation  (35) 

Hie  results  confirm  the  remarks  made  in  Section  II  about  the  relative 
worth  of  the  first  three  approximations,  i.e.}  only  (13)  has  merit 
unless  c = cry /a^  is  ve ry  close  to  .-no.  The  results  also  show  that 
(35)  is  closer  fco  the  exact  than  (13)  but  not  by  much. 

Let  us  now  consider  another  example  with  a smaller  c.  Consider 
Cy  » 15  and  r*y  = 100  so  that  e = .15.  In  this  case  y **  1.0450  with 
corresponding  multiplying  factor  .9714  obtained  by  interpolation. 


* 13  ft., 


100  ft. 


.6941  7jf 


Exact  ECEP  * .6941  ? « 69.1600  ft.,  Reference  13} 

Approx,  ECEP  - 1.1 774 08. 72m  - 45,6005  ft..  Equation  (12) 

Approx,  ECEP  » 1. 1774 (37. 5000)  * 67.7005  ft..  Equation  (13) 

Approx.  ECEP  « 1.1774(71. 5017)  * 84.1862  ft.,  Equation  (14) 

Approx.  ECEP  = .9714(71.3017)  * 69.4568  ft..  Equation  (35)  . 

Here  again  wo  e equation  (13)  is  closer  to  the  exact  than 
either  (12)  or  (14)  hut  barely  close  enough  for  practical  use.  On 
the  other  hand,  approximtien  (35)  is  quite  close  to  the  exact, 
differing  only  in  the  third  digit. 

Thus  tar  Che  eewent*  in  this  section  have  been  confined  to  the 
bivariate  case.  At  a for  emotioned,  the  need  for  such  an  approx  (mat  ion 
is  not  great  for  the  bivariate  case  since  we  have  tables  of  the  exact 
solution  in  Reference  [Si.  however,  when  the  number  of  d intensions 
exceeds  two,  tables  of  the  exact  solution  and/or  computer  programs 


for  the  exact  solution  are  not  believed  t~  exist.  Hence,  in  the 
general  multivariate  ease,  an  approximation  for  the  radius  of  the 
1Q0P%  k-sphere  could  be  of  practical  value. 

To  obtain  such  an  approximation,  we  apply  the  above  limiting 
concept  to  the  general  express  ion  for  L in  (34)  vice  the  two 
dimensional  expression  in  (24).  This  provides  us  with  the  following 
approximation  to  Sp,  the  radius  of  "he  IQQPX  k-sphere  in  the 
multivariate  case: 


t 

f (<7j  + ...  + o^)/kJ 


v 


+ . . . 


(36) 


The  worth  of  this  approximation  is  difficult  to  ascertain  due  to 
the  absence  of  tables  of  or  programs  for  the  exact  solution.  However, 
some  rough  comparisons  can  be  made  if  wo  once  mere  resort  to  simulation. 

To  make  these  comparisons,  we  specify  b,  and  tfien 

compute  our  approximate ion  to  ftp  wit ich  we  shall  designate  Sp.  For 
the  same  j,  Wi'  then  repeatedly  sample  front  a multivariate  normal 
popuUt ion,  each  time  sensing  or  the  radial  error  R.  The  estimate 
of'  P,  P,  is  simply  the  ratio  of  the  number  of  replicates  in  which  R 
Is  less  than  S,'  to  the  total  number  ef  replicates.  If  the^apprexi«a» 
tien  has  merit,  P should  lie  within  sampling  variation  of  p. 

Several  example  eases  are  considered  below,  the  reader  who  is 
seriously  interested  in  using  appro# last ion  (36)  should  engage  in  a 
•sere  extensive  evaluation.  Her  each  ease,  tire  number  of  replicates 
was  sat  at  10,000.  Using  this  number,  one  is  9$i  confident  that  the 
true  population  proportion  encompassed  b£  a U- sphere  of  radius  Sp 
ties  within  .01  of  s’.  Honte.  values  of  P within  .01  of  I*  Indicate 
that  the  approximation  Has  saerlt.  The  cases  considered  twt  the 
results  are  shown  in  Table  4. 
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TABLE  4 


Case  k P n, 

1 3 .50  l l 

2 3 .90  i 1 


3 3 .50  l 

4 3 .90  l 

5 4 .50  1 

6 4 .90  l 


r_2 

% 

S * 

P 

i 

1.5382 

.5031 

i 

2.54303 

.9033 

4 

3.6367 

.5032 

4 

7.1374 

.9020 

4 

3 

7.3596 

.4965 

4 

3 

14.3186 

.9019 

in  the  first  two  cases  in  Table  4,  die  distribution  is  spherical 
so  the  values  of  Sf\  for  chose ^eascs  are  exact.  We  note  that  P lies 
within  sampling  variation  of  P for  both  values  of  P ax  expected,  in 
the  last  four  cases,  the  distribution  is  ellipsoidal  and  as  such, 
values  of  3‘,  are  approx inat ions.  In  all  four  cases,  P lies  within 
sampling  variation  of  P,  because  of  the  email  number  of  eases 
examined,  we  cannot  conclusive Is-  state  that  the  approximation  is 
"good”  for  ail  eases.  However,  these  eases  are  evidence  that  the 
approximation  has  promise. 


VII.  CONCLUDING  REMARKS 
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The  addition  of  this  report  to  References  [7]  and  [8]  provides 
one  with  the  methodology  to  construct  (P,y)  tolerance  k-spheres  for 
the  uncorrelated  multivariate  normal  distribution.  The  spheres 
can  be  constructed  whether  the  sample  size  is  finite  or  infinite 
(corresponds  to  the  case  where  the  parameters  are  known)  and  whether 
the  standard  deviations  are  the  same  for  each  direction  or  different. 
References  [7]  and  [8]  address  the  case  where  the  standard  deviations 
are  the  same  for  all  directions  and  provide  the  exact  formula  for 
the  radius  of  the  k-sphere.  This  report  addresses  the  case  where 
they  are  different  and  provides  an  approximation  for  the  radius  of 
the  k-sphere.  This  approximation  provides  little  loss  in  accuracy 
for  the  cases  examined  even  for  small  sample  sizes.  Collectively, 
these  reports  provide  the  weapons  analyst  with  the  tools  to  make  a 
more  meaningful  assessment  of  the  delivery  accuracy  of  a weapon  than 
has  been  possible  in  the  past.  However,  one  closing  cautionary  point 
should  be  mentioned.  These  reports  assume  that  there  is  no  bias  in 
the  weapon  system,  i.e.,  that  the  long  run  average  or  expected  miss 
distance  is  zero  in  all  directions.  If  this  assumption  of  zero  bias 
cannot  be  made,  the  procedures  set  forth  in  this  report  and  in 
References  [7]  and  [8]  are  not  directly  applicable. 
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